Closed Geodesies and Billiards on Quadrics related to 
elliptic KdV solutions * 



We consider algebraic geometrical properties of the integrable billiard on 
a quadric Q with elastic impacts along another quadric confocal to Q. These 
properties are in sharp contrast with those of the ellipsoidal Birkhoff billiards in 
M™. Namely, generic complex invariant manifolds are not Abelian varieties, and 
the billiard map is no more algebraic. A Poncelet-like theorem for such system 
is known. We give explicit sufficient conditions both for closed geodesies and 
periodic billiard orbits on Q and discuss their relation with the elliptic KdV 
solutions and elliptic Calogero system. 

1 Introduction 

One of the best known discrete integrable systems is the billiard inside an n- 
dimensional ellipsoid (more generally, a quadric) 
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with elastic refections on The billiard inside Q (Birkhoff case) inherits the 

remarkable property of geodesies on Q given by the Chasles theorem: the trajectories 
or their continuations before and after impacts are tangent to the same set of n 
quadrics that are confocal to Q (see, e.g., jH]). The parameters of these quadrics 
deliver n independent and commuting first integrals of the discrete system. 

Veselov [3BJ described the Birkhoff billiard map in terms of discrete Lagrangian 
formalism and showed that its complex invariant manifolds are open subsets of 
coverings of Jacobians of hyperelliptic curves. Moreover, the restriction of the map 
to such manifolds is represented by shifts by a constant vector, which gives rise to 
an addition law on Abelian varieties. 

Explicit theta function solutions for this billiard were obtained in [HE] by applying 
spectral theory of difference operators, and later in ^S] by finding a degeneration 
of theta function solutions for geodesies on an ellipsoid of a higher dimension. 

The periodicity problem for the Birkhoff billiard map is closely related to a 
generalization of the following geometrical problem: given two smooth conies in 
the projective plane, C and D, to construct a closed polygon inscribed in C and 
circumscribed about D. The full Poncelet theorem states that, if a closed polygon 
exists for given C and D, then every point of C is a vertex of a closed polygon and 
all such polygons have the same number of sides |22| . 

Cayley theorem gives an explicit condition to find such a polygon. Namely, 
the iterative process of constructing the polygonal line is equivalent to shifts by a 
constant vector v in the Jacobian of an elliptic curve, £. Then the polygonal line is 
closed after n steps if and only if v is neutral on £ f |§1 12*3]). 

The relation of the Poncelet problem to periodic billiard orbits inside Q is evi- 
dent, and such orbits are described by generalizations of the theorems of Poncelet 
and of Cayley in QJ CHI H3 EH ■ 

On the other hand, one can consider geodesic motion on the ellipsoid Q combined 
with elastic reflections along its intersections with the confocal quadric 



d being an arbitrary fixed parameter. Thus one obtains the billiard on the ellipsoid 
described by the map B : (x,v) i— ► (x,v), where x,v £ M n+1 are respectively the 
coordinates of an impact point on Q n Qd and the outgoing velocity at this point, 
whereas x, v denote the same objects at the next impact point. 

As shown in ^3], this billiard map is also integrable and the trajectories between 
impacts (arches of geodesic) or their continuations are tangent to the same set of 
n— 1 quadrics (caustics) in M n+1 that are confocal to Q and Qd- For n = 2, algebraic 
geometrical features of this problem have already been discussed in [37] . 

Next, following ^3]; if a billiard trajectory of Q is closed, then all the trajec- 
tories sharing the same caustics are closed as well and have the same period and 
length. Thus one arrives at a generalization of the celebrated Poncelet theorem on 
a quadratic surface. 
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Contents of the paper. In this paper paper we study algebraic geometrical 
properties and periodic orbits of the billiard on the ellipsoid Q with impacts along 
its intersection with a confocal quadric Qd- 

All our statements in the complex setting also hold in the case of a general 
quadric Q. On the other hand, in the description of real closed geodesies and periodic 
billiard orbits we always assume Q to be an ellipsoid, since, for a generic real quadric 
Q, the billiard problem becomes meaningful only after a compactification. 

We show that the algebraic-geometrical properties of the billiards on Q with 
impacts along Qd are in sharp contrast with those of the Birkhoff case. Namely, the 
generic complex invariant manifolds are not Abelian varieties, but open subsets of 
theta-divisors of n-dimensional hyperelliptic Jacobians. 

Next, although the billiard map B is integrable and all of its first integrals are 
algebraic, in the complex coordinates (x, v) it is infinitely many valued, whereas 
in the real domain the new values (x, v) are determined by solving a system of 
transcendental equations which involve the inversion of hyperelliptic integrals. 

To our knowledge, such billiards actually provide a first nontrivial example of 
a discrete integrable system which cannot be described in terms of an algebraic 
addition law on Abelian varieties. 

Such systems can be regarded as discrete analogs of hyperelliptically separable 
systems, a class of finite-dimensional integrable differential equations whose generic 
complex solutions have movable algebraic branch points and are single-valued on an 
infinitely sheeted ramified covering of the complex plane t (see, e.g., |38| IH l2| IT8] ) . 

In Section 2 we briefly recall an algebraic geometrical description and theta- 
function solution of the Jacobi problem on geodesies on a quadric Q, then describe 
the main properties of the billiard on Q. 

In Section 3, we consider the periodicity problem for the billiard map B and we 
show that, in general, such conditions cannot be formulated in an algebraic form. 

This however becomes possible for special initial conditions for which the cor- 
responding hyperelliptic curve T is a covering of an elliptic curve £ and one of the 
holomorphic differential on V reduces to that on £. Under these conditions, the cor- 
responding geodesies on Q have a finite number of intersections with the boundary 
Q C\Qd and therefore are closed. 

In Section 4, we use known properties of the periodic KdV solutions and of the 
solutions to the elliptic Calogero system[I3 E3 El EEl 021 [HBl EI3 EH] to describe the 
curves T and the related closed geodesies on Q. As a result we show that to each 
elliptic A-soliton KdV solution or solution of the A-body elliptic Calogero systems 
satisfying the locus condition there corresponds a family of closed geodesies on a 
quadric Q. 

We consider in detail the case of the 3:1 coverings T i— ► £ and obtain explicit 
sufficient conditions for a real geodesic on a 2-dimensional ellipsoid Q to be closed. 

Finally, in Section 5 we derive sufficient conditions for the periodicity of the 
billiard map on Q by describing admissible boundary parameters d of Qd- To our 
knowledge, this is the first time a sufficient condition of periodicity for such billiard 
map is given. 
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In the conclusion we outline possible extensions of this approach to periodic 
billiard orbits on a multi-dimensional ellipsoid Q. 



2 Geodesies and Billiards on a Quadric 

We start with the celebrated Jacobi problem of the geodesic motion on an n- 
dimensional ellipsoid Q. The problem is well known to be integrable and to be 
linearized on a covering of the Jacobian of a genus n hyperelliptic curve. Namely, 
let I be the natural parameter of the geodesic and Ai, . . . , X n be the ellipsoidal co- 
ordinates on Q defined by the formulas 



* = /*- Al) ;-' (0 SH i = l,..,n + l. (2.1) 

Then, denoting X k = dX^/dl the corresponding velocities, the total energy ^(X,X) 
takes the Stackel form 

n 1 n 

H = 2 J2^7 X — ) U^- X j)^h *(A) = (A-a 1 )-.-(A-o n+1 ). 

k=l ^ k) j+k 

According to the Stackel theorem, the system is Liouville integrable and, after re- 
par ametrization 

dl = X 1 ---X n ds, (2.2) 
the evolution of Aj is described by quadratures 

A^ _1 dAi | | A*- 1 ^ = f ds for k = 1 , 



2V-*(Ai)i*(Ai) 2 y /-*(\ n )R(\ n ) [0 for k = 2,..., n 

R(X) = A(A — ci) ■ ■ ■ (A — Cn-i), 

Cfc being constants of motion. Then the components of the velocity have the form 

_ dXj _ ^Jai{ai- Xj) ■ ■ • (a— A n ) 1 y/ -V (X k ) R(X k ) 

The quadratures involve n independent holomorphic differentials on the genus 
n hyperelliptic curve T = {/j? = — \I/(A)-R(A)} and give rise to the Abel-Jacobi map 
of the n-th symmetric product r( n ) to the Jacobian variety of T, 



Pi Pn 

^>k-\ h J u k 

Po Po 



J Lo k = u k , k = l,...,n, (2.5) 



X k—1 dX / \ 



4 



where m,...,u n are coordinates on the universal covering of Jac(r) and Pq is a 
fixed basepoint, which we choose to be the infinity point, oo, on T. 

Since u\ = s+const and U2,---,u n =const, the geodesic motion in the new 
parametrization is linearized on the Jacobian variety of T. 

For the classical Jacobi problem (n=2), its complete theta- functional solution 
was presented in |39| . and, for arbitrary dimensions, in 25 , whereas a complete 
classification of real geodesies was made in [3]. 

Namely, let us fix a canonical homology basis on T and let ui\ , . . . , u n be the 
conjugated basis of normalized holomorphic differentials, 



n 



^ Cik Uk , (2.6) 



fc=l 



Cik being normalization constants so that the 2n x n period matrix of T has the 
form (27T/I, B), where B is the n x n Riemann matrix. 

Let tp = (ipi, . . . , ip n ) T , <fi = S"=i Ci S u s be the corresponding local coordinates 
on Jac(r) and 



(ip\B) = exp{{Ba, a)/2 + (<p + 2irj(3, a)} 9(<p + 2irjf3 + Bat), (2.7) 



be associated theta- function with characteristics a = (a±, . . . , a n ), (3 = (f3±, . . . , (3 n ) 
and for K,M G Z n , 



(ip + 2ttjK + BM) = exp(2vrj e) exp{- (BM, M) /2 - (M, ip)}e 

e = (a,K)-(P,M). 



(<P), (2-8) 



Then the inversion of the map (|2.5[) applied to formulas (|2.1j) leads to the following 
parametrization of a generic geodesic 

8[A + ri t ](U lS + ip ) 

Ms) = *i — , t = l,...,n+l, 2.9 

9[A\(Uis + ifo) 



U — 2{C\\, ... ,C n 



Here A = (5", 5'), rji = (r?f,r/^) G M 2n /21R 2ra are half-integer theta-characteristics 
such that modulo the period lattice of T the following relations hold 

K<h,o) 

torjrtf + Br£= / (wi, . . . , Cu n f, 2ttjS" + Btif = JC, (2.10) 

J CO 

/C being the vector of the Riemann constants, and K{ are constant factors depending 
on the moduli of T only. 

Remark 2.1. As follows from the quasiperiodic law (|2.8|) . expressions (|2,9|) are 
meromorphic functions on a covering Jac(r) of Jac(r) obtained by doubling some 
of its periods (see also [2*5]). 
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The billiard on Q. We now describe the billiard map associated to the geodesic 
motion of a point on Q with elastic bounces along the confocal quadric Qd- By 
Q2.1|) . when (Xi, . . . , X n+ \) £ Q Pi Qd, one of the points Pi = (Aj, /Xj) on T (without 
loss of generality we choose it to be P n ) coincides with one of the points Ed± = 
(d, ±y / R(d)). Under the Abel-Jacobi map (|2.5|) . the condition P n = Ed± defines 
two codimension one subvarieties in Jac(r), the theta-divisors 

e d ± = {9[A}(<pTCq/2)} G Jac(r), 

where C = {Co-) is the matrix introduced in (|2,6|) and q is the vector 

r E <i+ 

q = ( qi ,...,q n ) T = / K..., W „feC" (2.11) 
JE d _ 

@d+ and Qd- have a common (n — 2)-dimensional subvariety. 

Hence, for fixed constants of motion, the coordinates x of impact points on Qf~)Qd 
with velocities v are described by degree n — 1 divisors {Pi, . . . ,P n __i}, that is, by 
a point ip on Qd- or Qd+- 

An algebraic geometric description of the motion between impacts and at the 
elastic bounce along Q n Qd is given by the following proposition. 

Proposition 2.1 The geodesic motion on Q between subsequent impact points x 
and x on Qf]Qd corresponds to the straight line uniform motion on Jac(r) between 
Qd+ and Qd- along the u\-direction. The point of intersection with Qd- gives the 
next coordinate x and the ingoing velocity v' . At the point x the reflection v' — > v 
(from the ingoing to the outgoing velocity) results in jumping back from Qd- to Qd+ 
by the shift vector q, which does not change x. Then the procedure iterates. 



The process is sketched in Figure 1. 




Figure 1 



Proof of Proposition Yl.W First, we want to show that any subsequent impact points 
and the velocities (x, v) and (x, v ') correspond to intersections of the flow with Qd+ 
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and Od- respectively. For this purpose, we introduce the new time parameter H 
such that 

d " (A 1 -d)---(A n -d) dS - 



Notice that in the real case the number yj — ^/(d)R(d) is real and positive. Hence, 
as the point on Q moves from x to x, the parameter H changes monotonically from 
— oo to oo. On the other hand, in view of the quadratures ([2.3)1 . one has 



"J (Aj - d)y / -^(Xj)R(Xj) 



and, therefore, 



r p i J-y(d)R(d)dX 

/ ¥ + const. (2.12) 

7p (X-d)^-^(X 1 )R(X ] ) 



The latter expression contains a meromorphic differential of the third kind on the 
curve T having the pair of simple poles at E± = (d, ±y/—ty(d)R(d)) with residue 
±1 respectively, hence the sum ([2. 12)) is an Abelian transcendent of the third kind. 
According to the theory of theta- functions (see e.g., or [7]), under the Abel- 
Jacobi map ()2.5|) the sum takes the form 

„ , 0[A](^ + Cq/2) „ 

E = lo § fl A n /o + Cs + const ' 
6[A]((p - Cq/2) 

where C is a normalizing factor, C = (Cjfc) is defined in ()2.6)1 and the vector q is de- 
fined in H2.11JI . It follows that for 3 = — oo and S = oo the point ip in Jac(r) belongs 
to Qd+ an d @d— respectively, which proves the first statement of the theorem. 

Next, we notice that for any ip 1 £ @_ and (p £ 0+ such that (p = ip' + Cq, the 
formula ()2.9)) gives one and the same Cartesian coordinates X{. Indeed, under the 
mapping 1)2.5( 1. y> and are represented by the degree n divisors 

{Pi = (Ai, . . . , P n -\ = (A„_i,/%_i), P_} 

and, respectively, {Pi, . . . , P n -i, which by (|2.1|) . lead to the same values of Aj. 
On the other hand, expressions ()2.4)l imply that for (p and <p' the velocities are 
different. 

As a result, the intersection point ip' £ ©_ corresponds to the pair (x,v'), 
whereas ip = ip' + Cq £ + gives rise to x, and the outgoing velocity v. rj 

It follows that the restriction of the billiard map B : (x, v) — * (x, v) onto its 
invariant tori can be regarded as a transformation B c : Qd+ ~^ ©d+- Then, as a 
corollary of Proposition 12.11 we obtain 

Theorem 2.2 The map B c is given by the shift of coordinates it2, ■ ■ ■ ,u n by q2, • • • , q n 

respectively, whereas the first coordinate u± on Jac(r) is chosen such that the vector 
u = B c {u) remains on the theta-divisor &d+- 

Indeed, according to Proposition 12 . 1 1 (we recall that ip = Cu), under the straight 
line flow on the Jacobian from Qd+ to Qd~ , the coordinates U2 , • • • , u n remain un- 
changed, and under the passage to 0^+ they increase by q2, • • • , q n - 
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Remark 2.2. For n = 2, the above algebraic geometrical property has already 
been described in [37]. The algebraic geometric description of the billiard on Q 
stands in sharp contrast to the one of the Birkhoff billiard system - as well as to 
most of the known integrable maps. Indeed in the Birkhoff case, the billiard map is 
given by a shift by a constant vector on a complex invariant manifold which is an 
open subset of an Abelian variety. In our case, the shift takes place on a non-Abelian 
variety and, ctS Si result, it is not fixed. 

Branching of the map B c . The above properties of the map B cannot be used to 
derive an explicit expression for the billiard map (x,v) — ► (x,v), since the complex 
map £> c : Qd+ —> ©d+ described by Theorem 12.21 is infinitely many- valued. Indeed, 
let ®d+ be the universal covering of 0^+ in C n = (ui, . . . , u n ) and let tt denote the 
natural projection @d+ — > C n_1 = (u2, . . . ,u n ). 

Theorem 2.3 (see also 1 ). Suppose that the curve T does not have nontrivial 
symmetries and its Jacobian does not possesses Abelian subvarieties, then 

1) . Under the projection 7r, @d+ is an infinitely sheeted covering o/C n_1 ramified 

along a codimension one subvariety of &d+ ■ 

2) . A complete preimage of any point in C n ~ l consists of an infinite number of 

points on @d+ > which, in turn, represent an infinite number of points on ®d+ ■ 

Proof. Let Oi, . . . , t>2n £ C n_1 = (u2, ■ ■ ■ , u n ) be the projections of 2n independent 
period vectors of Jac(r). For any point V 6 @d+, the projections of points of its 
equivalence class on have the form 



Let z be any fixed point on C™ . Under the condition of the theorem, the projec- 
tions t)j are incommensurable and the integer coefficients rrij may always be chosen 
in such a way that the above sum represents a point in any small and a priori de- 
fined neighborhood of z. In this sense, any point on + is projected into such a 
neighborhood of any point on C n_1 . Hence, a complete preimage 7r _1 of any point 
on C n_1 corresponds to an infinite number of points on the theta-divisor. 

In the particular case n = 2, when the theta-divisor or its translate coincides 
with the genus 2 curve T itself, the above property is known in connection with the 
inversion of a single hyperelliptic integral 



Namely, as was shown by Jacobi (see, e.g., }28J ) , the covering r — » C = {ui} has 
infinitely many branch points whose projections onto the complex plane u\ form a 
dense set. 




2n 
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Hence, given fixed coordinates U2, ■ ■ ■ , u n on Jac(T), there is an infinite number of 
Ui, which define the next impact point x S QdQd via Proposition l2.il In particular, 
this implies that the billiard map B : (x, v) — ► (x, cannot be algebraic. 

The above phenomenon is also explained by the fact that a generic geodesic on 
Q intersects any connected component of the boundary Q n Qd at infinitely many 
points. In the real domain, given x G Q Pi Qd) it is possible to determine x by 
choosing the nexi intersection with the boundary, but in the complex setting such 
a choice cannot be made. 

Remark 2.3. Due to the above properties, the billiard map B can be regarded 
as a discrete analog of the hyperelliptically separable systems studied, in particular, 
in [HHIEEIE])- Such systems, although being Liouville integrable, are not alge- 
braic completely integrable: their generic complex solutions have movable algebraic 
branch points and are single- valued on an infinitely sheeted ramified covering of the 
complex plane t, so they possess the so called weak Painleve property. The branch- 
ing of the billiard map B can be viewed as a discrete version of such a property (see 
also |2H). 

3 Periodic Billiard Orbits on Q D Qd 

We now concentrate on algebraic conditions for a billiard trajectory on the quadratic 
surface Q to be periodic. 

As mentioned in the Introduction, the Birkhoff billiard in the ellipsoid Q C 
is described by translations on the Jacobians of genus n hyperelliptic curves 



Then the corresponding billiard trajectory is periodic if, for a certain integer m, the 



J Bo- 
is neutral in Jac(r). The periodicity conditions written explicitly in terms of the 
coefficients of the curve T were obtained in ^j] by proving the following theorem. 



r = {y 2 = -R(x)} Cl 2 = (x, y), R{x) = (x - h) ■ ■ ■ (x - b 2n+1 ). 



vector 



m 




Theorem 3.1 For a constant d £ C and involutive points Ed± 

r, the vector m j^ d+ oj (m > n) is neutral in Jac(r) if and only if 





where the entries Sj are determined from the expansion 



\/(x - 61) • • • (x - b 2n +i) = Sq + Si(x - d) + S 2 {x — d) 2 H . 
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This theorem provides a multidimensional generalization of Cayley's algebraic con- 
dition for the existence of a polygon inscribed in a conic and circumscribed about 
another conic. 

For the billiard on the quadric Q, the situation is quite different: generic com- 
plex invariant manifolds are no more Abelian varieties, but theta-divisors. As fol- 
lows from Theorem 12.21 the billiard trajectory on Q is periodic if for some integers 
m, mi, . . . , rri2n the following equality holds 



Since the projections are generally incommensurable, the above condition can always 
be fulfilled, hence in the complex setting this criterion for closed billiard trajectories 
makes no sense. 

Remark 3.1. On the other hand, it seems natural to try to replace the above 
condition with the following relation in the Jacobian: 



That is, although the map B c is not given by a shift by q on Jac(r), it is required that 
after N iterations it yields an identity map on Jac(T) and therefore, on the theta- 
divisor 0^+ as well. This condition was proposed, amongst other results, in the 
recent paper [I^], and explicit algebraic conditions on the coefficients of the curve T 
were derived from it. However, since a generic geodesic X(s) C Q intersects QCiQd 
at infinitely many points, these conditions do not guarantee that the corresponding 
closed billiard trajectory consists of pieces of geodesies joining the nearest points of 
intersection with Q n Qd- 

Periodic billiard orbits related to closed geodesies. The projections X>i, . . . , X>2n 
may become commensurable and the closeness condition in the complex domain 
may become meaningful if the hyperelliptic curve T has certain symmetries and its 
Jacobian contains Abelian subvarieties, as described by the Weierstrass-Poincare 
theory of reduction of Abelian varieties and functions (see e.g., [27] and, for modern 
exposition, [HE])- 

In this paper we concentrate on the case when the genus n hyperelliptic curve T n 
is an A r -fold covering of an elliptic curve £ and the "first" holomorphic differential 
u\ = d\j \x reduces to a holomorphic differential on E. 

Under these conditions, the ui-flow on Jac(T) intersects @d+ at a finite number 
of points. As a result, the corresponding geodesies on Q also have a finite number 
of intersections with the boundary Q n Qd and therefore are closed. 

According to the Poincare theorem, the existence of covering r h f implies 
certain conditions on the components of the Riemann matrix of the curve T (0 
E]). However, our objective below is to give sufficient conditions in terms of the 
coefficients of the polynomial R(X). 



2;i 
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Remark 3.2. We emphasize that not any covering r i— > £ corresponds to closed 
geodesies. For example, in the simplest possible case n = 2, N = 2 studied by Jacobi 
and widely described in the literature (see, e.g. H3|7j), the curve V is birationally 
equivalent to the canonical curve 

w 2 = z(z — l)(z — a)(z — /3)(z — a/3), 

a, (3 being arbitrary constant. It covers two elliptic curves 

l ± = {Wi = Z ± (l - Z ± )(l - k ± Z ± )}, kl - T ^ ? 



(l-a)(l-/?) 



and both of the holomorphic differentials u\ ,uJ2 on T reduce to linear combinations 
of the holomorphic differentials on £ + and t— 

In this case the inversion of the quadratures ()2.3|) and formulas (|2.9j) lead to 
solutions in terms of elliptic functions of £±, but, since their periods are generally 
incommensurable, the corresponding geodesies are quasi-periodic. Such geodesies 
were already studied in [Sj. 

4 Elliptic KdV solutions, elliptic Calogero systems and 
closed geodesies on a quadric 

The most classical series of hyperelliptic tangential coverings Q n — > £ (see [21]) arises 
in the spectral theory of the Lame potentials or as a subclass of spectral curves of 
the elliptic Calogero (EC) systems, which describe the motion of N particles with 
coordinates qi,...,qw on the circle under the potential J2i<j p(Qi ~ Qj)i see CHI- 
Here p{q) = p(q | uj,uj') denotes the Weierstrass elliptic function with half-periods 
iv, (J . 

The EC systems are integrable and, as shown in they possess a matrix 
N x N Lax pair with an elliptic spectral parameter u G £ = C/{2lu7j + 2u/Z} 1 , 

j t L(u) = [L(u),M(u)}. 

The spectral curve Q = \L(u) — fjl\ = (the so called Krichever curve) is an iV-fold 
covering of the torus £ and its genus does not exceed N. 

The EC systems possess invariant subvarieties given by the KdV locus condition 

N 

^2p'(qi-qj)=0, qi^qj, j = l,...,N. (4.1) 

If N equals a "triangular number", i.e., N = n(n + l)/2 for some n, and the 
coordinates qi satisfy the KdV locus condition, then the following properties hold 

We do not give the explicit expressions for L(u) here 
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1). ( [171 E|) For a certain constant C, the function 



N 



U(x,t) = 2^p(x- qi (t)) + C (4.2) 



i=l 

is an elliptic solution of the KdV equation Ut = 6UU X — U xxx . 

2) . ([32, 33 ) The Krichever spectral curve is hyperelliptic of genus n: there exists 

a birational change of variables (p(u),n) — > (z,w) such that the equation 
defining Q takes the canonical form 

W 2 = R2n+l(z) = (Z — Zl)---(Z- Z2n+l) (4.3) 

with some parameters Z\~. Up to a constant factor, the last holomorphic dif- 
ferential oj n = z n ~ l dz/w on Q reduces to the holomorphic differential dp/p' 
on £. 

3) . ( 26 ) Let ipi, . . . , i/j n _i,i/j n be the local coordinates on Jac(^) associated to the 

holomorphic differentials u>k = z k ~ 1 dz/w, k = 1, . . . , n on Q and let (pi, . . . , (p n 
be the coordinates associated to the normalized holomorphic differentials. Let 
9 ((f) be the theta- function associated to Q and U G C n be the tangent vector 
of the curve Q C Jac(C?) at oo G G, which in the coordinates Vfc has the 
form U = (0, . . . , 0, 2) T . Then for suitable constant vectors V, W £ C n , the 
transcendental equation 

6[A](Ux + Vt + W) = 0, q,teC, (4.4) 

has exactly N solutions x = q%(t), . . . , x = quit) on Jac(C7), which provide the 
solution for the EC system. 

4) . ( 20 ) One can choose a canonical homology basis on Q and the basis of nor- 

malized holomorphic differentials such that the vectors U, V in (|4.4[) take the 
form 

U = (1,...,0,0) T , V=(0,V 1 ,...,V n ^ 1 ) T (4.5) 
and the period matrix of Q becomes 



(2irjIB), B 



( t/N 2ttj/N ••• 0\ 
2-kj/N 


V o / 



(4.6) 



where (2irj, r) is the normalized period matrix of the elliptic curve £ and 23 
is a (n — 1) X [n — 1) Riemann matrix. 

In view of normalization (|4.5|) . the function 9[A](Ux + W) can be written as 
9[A](x,t), t = (t2, ■ ■ ■ ,t n ). It is a one-dimensional theta-function of order N 
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with respect to x and it admits the factorization 

N 

9[A](x,t) = C(t) Y[O n (x - «j(t)), (4.7) 



where 9u(x \ t) = 



1/2 

1/2 



3=1 



t) is the one-dimensional Riemann theta- 



function associated to £ and the factor £(t) depends on t only. 

Notice that item (3) provides an implicit solution for the EC system. Explicit 
expressions for certain symmetric functions of the coordinates qt were first given in 

noi. 

According to the Poincare reducibility theorem (see e.g., [7J, apart from the 
curve £, the Jacobian of Q contains an (n — l)-dimensional Abelian subvariety A n -\ 
and is isogenous to the direct product £ x A n -\. 

Notice that for n = 2, explicit algebraic expressions of the covers and coefficients 
of hyperelliptic curves are known for N < 8 (see [35] ) ■ 

Remark 4.1. Property 3) also implies that the x-flow (or ip n -&ow) on Jac(<5), 
which is tangent to Q C ,lac(G) at its infinity point oo, intersects any translate of the 
theta-divisor O C Jac(£?) precisely at n(n + l)/2 points (possibly with multiplicity). 



This property is similar to what we require for the linearization of the geodesic 
flow on the re-dimensional quadric Q. However, as seen from the quadratures ()2.3j) . 
the linearized flow on Jac(r) has a different direction: it is tangent to the embedded 
hyperelliptic curve T C Jac(r) not at oo, but at the Weierstrass point O = (0,0). 
Nevertheless, both behaviors are equivalent: if (/?, 0) is a finite Weierstrass point on 
Q, then by a birational transformation 

. a w 

X = J^WY ^ = (z- ( } 

the points oo and (/3, 0) on Q are sent to O and oo on T, respectively. Then, 
identifying the curves Q and T as well as their Jacobians, we find that the KdV flow 
(V'n-Aow) on Jac(<5) is represented as ui-flow on Jac(r) and vice versa. 

Thus, the birational transformation ()4.8|) establishes a relation between the el- 
liptic iV-soliton KdV solutions (|4.2|) or the solutions of the A r -body elliptic Calogero 
systems satisfying the locus condition and the geodesies on a quadric Q that are lin- 
earized on Jac(r). Moreover, since the KdV solutions are periodic, all such geodesies 
are closed. 

Clearly, for the geodesic and the quadric (ellipsoid) to be real, certain reality 
conditions on the coefficients of Q must be satisfied. 

In particular, assume that the roots Zi, . . . , ^2n+i of Q (see formula (|4.3[l ) have 
such values that the above transformation with certain a and (3 = Z2 n +i sends them 
to real positive numbers ei, . . . , e2 n and the infinity. In this case one can always 
choose a partition 

{e%,. . . ,e 2n } = {ax, ■ ■ • ,«n+i} u { c i> • • • ,Cn-i} 
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such that the quadratures (|2.3j) describe the real closed geodesies on the ellipsoid Q 
that are tangent to the confocal quadrics Q(ci), . . . , Q(c n -\). 

Parametrization of the closed geodesies. Let the hyperelliptic curve T be 
isomorphic to Q via (|4.8|) . Then, in view of property (4) of the EC solutions, we 
choose local coordinates (p±, . . . ,tp n on Jac(r) such that its period matrix has the 
form (|4~6|) and the function 0[A](Uxs + cpo \ B) defined in (|2.9|) . (|2.10j) can be written 
as #[A](s,t | B), t = (t2, ■ ■ ■ ,£ n ) T =const. Here, as above, s is the new parameter 
introduced in (|2.2|) and t plays the role of the constant phase vector numerating the 
geodesic. 
Now let 

{qi(t), . . . ,qN(t)} and {pii(t), ■ ■ ■ ,pm(t)}, i = l,...,n 

be solutions to equations 

0[A]O,t) = O, respectively 0[A + rji](s, t) = 0, seC. (4.9) 

Proposition 4.1 1). The Cartesian coordinates Xj of the closed geodesies admit 
the following parametrization in terms of theta- functions of 8: 

t) = ft t exp (s/2) • — — — j-—. (4.10) 

i/ere the factors £«(t) depend on the phases only, and, for any t, modulo the 
periods (2ttj,t) of 8, 

Pu + ---+Pm = qi + --- + qN + r/2. (4.11) 

2). Let the parameter d ^ {oi, • • • , a n +i, c±, . . . , c n _i} 6e such that the real closed 
geodesies on Q have a non-empty intersection with the quadric Qd- Then the 
geodesies intersect each connected component Ddj of Qn Qd at least at 2 and 
at most at n(n + 1) distinct points. 

As follows from (|4.1()|) . (f4. 1 1|> . as well as the quasi-periodic law of #n(s), when 
the argument s changes by the period 2itj of £, the coordinates X{ change sign, and 
they remain unchanged under the shift s — > s + r. Hence the corresponding closed 
geodesic is a 2-fold unramified covering of £ , i.e., it is an algebraic curve in M n+1 . 

Proof of Proposition 14.11 1). Due to definition of generic theta- functions with 
characteristics (12.7(1 . one has 

9[A + th](<p) = exp{(ri Brfo + (<p + 2^(r?f + A"),^)} 0[A](<p + 2njrf! + Brfo 

and therefore, 



e[A + m\(s,t) = Xi(t) exp(s(7 ?i )i) 

■9[A] L + r^ + 2^ 



N 



+ {m)'i 



,t + t (4.12) 
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with certain factors Xj(t) and vector t G C™" 1 . One can show that, for this choice of 
the basis of cycles and of the normalized differentials on T, one always has (f/iX = 
1/2. Next, in view of (|4.9|) and (|4.7|) . the theta- function in the right hand side 
of (|4.12|) admits factorization in one-dimensional theta- functions 6u(s — Pki \ t). 
Now applying this to the solution (|2.9|) for the Cartesian coordinates of a generic 
geodesic, one arrives at expressions (j4.1U|) . 

Since the geodesies are closed, these expressions must be doubly periodic in s. 
Then, due to the quasiperiodic law for 9n(s), relation (|4.11|) must hold. 

2). As shown in Section 2, under the Abel-Jacobi map (|2.5|) the intersections of 
the geodesic with Q n Qd correspond to points on the two translates @d~,@d+ of 
the theta-divisor on Jac(T). Next, according to property 3) of the EC solutions, the 
til-flow on Jac(r) intersects any translate of O at N = n{n + l)/2 (complex) points. 

According to Remark 2.1., the geodesic flow is linearized on a covering Jac(T) 
obtained by doubling of some of its periods. As a result, the ui-flow on Jac(T) has 
in total 4iV = 2n(n + l) intersections with the corresponding coverings of Od-L)Qd+, 
which gives rise to the same number of complex intersections of the geodesic with 



For d aj, the real manifold Q n Qd consists of 2 connected components 
Ddi,Dd2- Then, due to symmetry, each of the components is intersected at most 



The case N = 3, n = 2. In this simplest case the connected components Di d, D<i,d 
are ovals which bound simply-connected domains T>i t d,T^2,d on the 2-dimensional 
ellipsoid Q. Then one can define ingoing and outgoing intersections of the geodesic 
with Di d (when the geodesic goes inside, respectively outside of In view of 

Proposition I2.1[ the former ones correspond to the intersections of the Ui-flow on 
Jac(r) with 6^+, whereas the latter correspond to the intersections with the other 
translate 8^-. 

As follows from item (2) of Theorem 14. 1| the closed geodesic has 3 ingoing 
and 3 outgoing complex intersections with each component Di d- If d — > a>i, both 
components are confluent to the ellipse Q n {Xi = 0} and, therefore, the closed 
geodesies intersect the hyperplane JT, = at 6 (generally complex) points. Since 
the geodesic is an elliptic curve in C 3 , it can be represented as an intersection of the 
quadric Q with a cubic surface ( 30 ). 

The tangential 3:1 covering Q i— ► £ is associated to 3-elliptic KdV solutions 
and dates back to works of Hermite and Halphen f|24|). The genus 2 curve Q is 
birationally equivalent to 



QnQ d . 



at n(n + 1) points. □ 



H = <w 2 = - 




(4.13) 



which covers the elliptic curve with moduli (72,53, 



£ 1 = {W 2 = 4Z 3 



52^-53} C (Z,W). 
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The covering is given by the relations 

1 z 3 - 27 g3 2 W (z 3 - 9g 2 z + 54g 3 ) 

9 z2 - 3 52 ' " 27 - 3 52 )2 ' 

and possesses 2 branch points. The last (second) holomorphic differential on H 
reduces to an elliptic one: 

z dz 2 dZ 



w 3W 

In this case there is another 3:1 covering IT2 : H 1— ► £ 2 of the elliptic curve 

£2 = {W 2 = Z 3 - G 2 Z - G 3 }, (4.14) 
27 243 

G , 2 = -(5 2 3 + 9 53 2 ), £3 = — ^ai-al), (4-15) 



16 wz ao " 32 



such that 



I,, . /o3\ dz 1 cii? 

Z = --(4z 3 - 9 52 z-9(73), W = -i« z 2 -- ff2 , and — = -— . (4.16) 



4 V " ou " V 4 / w 3 W 

Let (iti,u 2 ) be local coordinates on Jac(-ff) associated to the holomorphic dif- 
ferentials (ki zdz/w,k2dz/w), k\,k2 =const. According to property (4) of the el- 
liptic KdV solutions, for suitable normalization constants k\,k2 the period matrix 
of Jac(if) takes the form 

2ir 3 r/3 2ttj/3\ 
2ttj 2ttj/3 b/3 J ' 1 j 

where (27rj, r) and (2-7TJ, b) are the normalized period matrices of £\ and £2 respec- 
tively. 

Algebraic conditions for the closed geodesies. Under the birational equiva- 
lence of the curves H and T, the structure of H imposes conditions on the parameters 
{a±, 0,2-1 03, ci} = {ei, . . . , e^} (squares of semi-axes and the constant of motion), for 
which the geodesies are closed, and which can be found explicitly. Notice that one 
can always set e± = a\ = 1 choosing appropriately a in (|4.8|) . Then 00 £ T is 
transformed either to one of roots ±^/3g2 or to one of the roots of Az 3 — 9g2Z — 53 
in (|4.13j) . In the first case e 2 , es, are positive, hence Q is an ellipsoid, in the other 
case some of are negative, so Q is a hyperboloid. In the latter case, reality condi- 
tion implies that is the geodesic flow on the hyperboloid corresponds to unbounded 
motions, which we do not consider here. 

These observations lead to the following sufficient conditions for a real geodesic 
on Q to be closed. 

Proposition 4.2 For given parameters e\ = 1, e.2 > the geodesic associated to 
the covering H — > £ is closed if and only if one of the following pairs of conditions 
on e3 and e^ is satisfied 

(l + i + J-Y 2 -4(i + i + ^-)=0, ! e2S l y MS) 

V e 2 e 3 y \ e 2 e 3 e 2 e 3 y 2(e 2 + e 3 + e2e 3 ) 
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or 

9e|e|-24e 2 e3(e 2 +e3) + 16e 2 e3 + 16e| + e| = 0, e 4 = ■ r. (4.19) 

3e 2 e 3 - 2(e 2 + e 3 ) 

If \4.18\) is satisfied, then the following ordering holds 

for < e 2 < -, < e4 < e 2 < e3 < 1, 
1 

for - < e 2 < 1 < e4 < e 2 < 1 < e3, 

for 1 < e 2 < 4, < e 4 < 1 < e 2 < e 3 , 
for e 2 > 4, < e4 < 1 < e3 < e 2 . 

Condition \4-19j) is satisfied if and only i/e 2 > 1, and the following ordering of roots 
holds: 

4 

for 1 < e 2 < -, < 1 < e 2 < e3 < e4 
4 

for - < e 2 < 4, < 1 < e% < e 2 < e 4 , 
for e 2 > 4, < 1 < e3 < e 4 < e 2 . 

Remark 4.2. The first condition corresponds to the case when e\ = 1 is trans- 
formed into one of the roots of the polynomial 4z 3 — 9g 2 z — 53, while in the second 
case 1 is transformed to y / 3g 2 . 

The inverse relations are described by the following lemma. 
Lemma 4.3 //e 2 ,e3,e4 is a solution of 14.18Q, then 



f3 = -Ul + ± + ±-), 92 = U\ 53 = -l( / 3 + l)^ + lV / 3+l 
3 V e 2 e 3 / 3 27 V e 2 / \ e 3 



If condition \4-iyty is satisfied, then 



1 1 1/11 

P = --, 92 = —, 93 = -r= - 



2 " 12 " 27 V8 e 2 e 3 e 4 ^ 
Proposition 14.21 and Lemma 14.31 are proved by pure calculations. 

Remark 4.3. Since Proposition 14.21 imposes two independent conditions on the 
parameters o 2 ,a3,c, not any triaxial ellipsoid has closed geodesies associated to 
the 3:1 tangential covering H — * £ . Naturally, this does not exclude existence of 
closed geodesies for the other types of hyperelliptic tangential coverings considered 
in Description of such geodesies will be the subject of a future study. 
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5 Periodic billiard trajectories on a 2-dimensional quadric 



As shown in Section 2, a criterium for a generic billiard on an ellipsoid Q to be 
periodic cannot be established in the complex setting. Below we restrict ourselves 
to 2-dimensional ellipsoid Q for which the trajectories between impacts on QnQd are 
segments of the closed geodesies X(s) described by Proposition 14.21 We study the 
following problem: given a family of such closed geodesies with an appropriate fixed 
caustic parameter c, for which boundary parameters d the corresponding billiard 
orbit is periodic? 

Correctness of the problem. As stated in the previous section, the closed 
geodesic X(s) has 3 ingoing and 3 outgoing complex intersections with each com- 
ponent D i c i of Q fl Qd- Hence, to an initial set (x,v), x £ there correspond 
3 possible branches (x,v), i.e., the restricted complex map B c is 3-valued. In this 
context the problem of its periodicity makes no sense, unless other properties are 
taken into account. 

In the real domain there may exist one, two, or three real branches {(x,v)}, 
and one naturally chooses a unique set (x*,v*) £ {(x,v)}, such that x and x* 
are subsequent points of intersection of X(s) with Di t d- However, (x*,v*) cannot 
be distinguished form the other branches by purely algebraic methods. Yet the 3 
branches have an important common property. 

Lemma 5.1 Given an initial set (x,v), the 3 sets (x,v) generate one and the same 
closed geodesic X(s) on Q. 

Indeed, according to Theorem l2.21 under the billiard map B c the local coordinate 
U2 on Jac(r) increases by the integral (\2 = Je^ + ~If> which depends on d only. 
Hence, the images of the sets (x, v) in Jac(T) have the same coordinate v>2 and, 
therefore, they belong to the same trajectory of the ui-flow on Jac(r). rj 

It follows that the map B c induces a map on the set of the closed geodesies, 
T c ^ : X(s) 1— > X(s), which is one-to-one. 

If the orbit of J- C) d is periodic, then the orbit of the transcendental real one-to- 
one map B* : (x,v) — > (x*,v*) is periodic as well. Hence, an algebraic condition of 
periodicity for T c ,d gives us a correctly defined periodicity condition for the billiard 
on Q. 

Periodicity conditions for J- C d- The map J- C d can be described as the shift 
U2 1— > U2 + c|2 and under the birational transformation i)4.8j) the increment q2 takes 
the form 
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Next, under the second covering II2 : H £ 2 C (Z,W) described by (|4.16[) . the 
above integral reduces to the elliptic one: q 2 = \l p , 

l p = / = , U 2 (E S± ) = (p, ±Vp*-G 2 p-G 3 ) , 

Jn 2 (E s _) yZ A - G 2 Z - G 3 

where 

p = -^-9/Ag 2 5-9/4g 3 ). (5.1) 

As follows from the structure of the period matrix of Q in (|4.17|) . the shift 
u 2 1— »• u 2 + c\2 results in the same ui-winding on Jac(r) if c\ 2 belongs to the lattice 
Ao = {27rjZ + ^bZ}. Hence, the orbit of J- c ^ is periodic if and only if q 2 is a finite 
order point of Ao. 

Since the period lattice A = {2-7T/Z + bZ} of Jac(£ 2 ) is a sublattice of Ao, we 
conclude that T c ^ is periodic if, in particular, for an integer m, the integral ml p is 
neutral in Jac(£ 2 )- (Notice that one should exclude trivial cases, whenZp is neutral 
in the lattice Ao-) 

Then, by applying the Cayley-type condition (Theorem I3.1j) . one obtains an 
algebraic equation on the parameter p. 

In practice, for large m, such equations appear to be rather tedious. We consider 
only the simplest nontrivial case g = l,m = 3, when the condition (|3,1|) for the 
elliptic curve £ 2 takes the form 



54 S3 

55 S4 

Sj being the coefficients of the expansion 



= 0, (5.2) 



V 7 ^ 3 - G 2 Z - G 3 = y/(Z - pf + P 2 (Z - pY + P l (Z-p) + P 
= Sq + S^Z - p) + S 2 (Z - pf + ••• , 

where 

Po = P 3 ~ G 2P - G 3 , Pi = 3p 2 -G 2 , P 2 = 3p. (5.3) 
Condition ()5.2I) yields the equation 

(4P P 2 - Pi) 

■ (512P 4 + 3Pf - 384P 3 PiP 2 - 20P PiP 2 + 96P$P? + 64P 3 P 2 3 + 16PjfP?P%) = 0. 

Substituting here expressions (|5.3|) . we conclude that the integral ml is neutral in 
Jac(£ 2 ) if p is a root of one of the following two equations 

3p 4 - 6G 2 p 2 - 12G 3/ 9 - G| = 0, (5.4) 
- p 12 + 22G 2 p 10 + 22OG3 p 9 + 165G| p 8 + 528G 2 G 3 p 7 
+ (1776G^ - 92G 3 2 ) p 6 + 26AG 2 2 G 3 p 5 + (185G| - 960G 2 G!) p 4 
+ (80G^G 3 - 320G 3 ;) p 3 + (62AG 2 2 G 2 3 - 90G 5 2 ) p 2 + (896G 2 G| - 132GtG3) p 
+ 512G^ + 3G 6 2 - 96G 3 2 G 2 3 = 0. (5.5) 
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Remark 5.1. The fact that equation 1)5.2(1 has 16 roots {pj} is natural: they 



correspond to 32 points j-E^i = [Pji^yPj ~ &2pj — G3J j on the curve £2 that 
satisfy conditions 



O being a Weierstrass point on £2. The corresponding 32 points A(E Pj ±) £ Jac(<?2) 
form a complete set of sixth order points {(Z/6Z)2-7rj + (Z/6Z)b/2} of which the 
four second order points {0, ir j, b/2, n j + b/2} are excluded. 

The resulting algorithm of constructing real periodic billiard trajectories on Q 
consists of 4 steps. 

1) Given a family of closed geodesies on Q defined by parameters e\ = 1, &21 63, £a 

in Proposition 3.2, one finds the constant (3 in (|4.8|) and moduli 52^53 of the 
first elliptic curve £\ by using Lemma 14.31 and then, in view of (|4.15|) . the 
parameters G2 , G3 of the curve £2 ■ 

2) One finds the set {pj} of all roots of equations (|5.4j) . (|5.5() and then rejects the 

"trivial" ones, for which the integral 2 p is neutral in the lattice Ao. 

3) One computes the sets of corresponding admissible parameters 5 and d by solving 

the cubic equation (|5.1I) and the linear equation 5 = 1/d + (3. 

4) One selects the subset of real positive admissible values of d for which the quadric 

Qd C M 3 has a nonempty intersection with the closed geodesies on Q, that is, 
the values that belong to segments {— A(A — e\) ■ ■ ■ (A — e^) > 0}. 

The parameters e\ = I,e2,e3,e4 and the selected real values of d characterize 
the billiard boundary Q n Qd and the caustics for which the nontrivial billiard orbits 
are periodic. 

Conclusion 

In this paper we studied algebraic geometrical properties of the billiard on a quadric 
Q with bounces along its intersection with the confocal quadric Qd- 

If Q is a 2-dimensional ellipsoid, we explicitly described a special class of closed 
geodesies associated to the 3-soliton solutions of the KdV equation corresponding 
to the 3- fold tangential covering Q 1— > £ 1. 

The existence of the second elliptic curve £2 C Jac(£?) in this case enabled 
us to apply the Cayley-type periodicity condition (Theorem 13. 1J) as an algebraic 
description of finite order points on £2 and, thereby, to obtain an explicit description 
of 3-periodic billiard orbits on Q(~)Qd whose points are joined by arches of the closed 
geodesies. As mentioned above, for period m > 3 the periodicity conditions lead to 
quite tedious algebraic equations for the boundary parameter. 

This approach cannot be directly extended to construct similar-type periodic 
billiard orbits on an n-dimensional quadric Q, since, instead of the complimentary 
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elliptic curve £2, the Jacobian of Q contains an (n — l)-dimensional Prym subvariety 
A n -i, which, in general, is not a Jacobian variety or its covering (except n = 3). 

Thus, we arrive at the problem of an algebraic description of finite-order points 
in A n -i in terms of points of the curve Q, which deserves a separate studying. 
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